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Abstract 

We calculate the isovector response function of hot nuclear matter using various effective Skyrme 
interactions. For Skyrme forces with a small effective mass the strength distribution is found to be 
nearly independent of temperature, and shows little collective effects. In contrast effective forces 
with an effective mass close to unity produce sizeable collective effects at zero temperature which 
disappear at temperatures of a few MeV. We discuss the relevance of these results for the saturation 
of the multiplicity of photons emitted by the giant dipole resonance in hot nuclei beyond T=3 MeV 
observed in recent experiments. 
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The purpose of the present paper is to investigate the interaction dependence of the response 
function of hot nuclear matter to a small isovector external interaction. The motivation for such a 
study is to clarify why some calculations of these resonances using Random Phase Approximation 
(RPA) predict a sharp peak at zero temperature and an increase of the width with temperature [ll |2| 
whereas such is not the case for other self consistent RPA calculations |J . In this last case one already 
has at zero temperature a broad resonance which shows little evolution with temperature. In reference 
0] it was found that the collectivity and temperature dependance of the nuclear matter response is 
quite sensitive to the strength of the particle-hole interaction. This result was however obtained in 
the case of a schematic Skyrme type interaction which contains no momentum dependent terms and 
therefore cannot be considered as a reliable effective force. For this reason it appears necessary to 
perform again such a calculation using a full Skyrme force. Although calculations for a complete force 
are somewhat cumbersome we will show that analytical and transparent formulae can still be obtained. 
These formulae generalize those of Garcia- Reccio et al. |j] for the zero temperature response function. 

To calculate the response function at finite temperature in nuclear matter, we consider an infinites- 
imal external field of the form 

V ext = er 3 e- i ^e- i ^ t , (1) 

where T3 is the third isospin Pauli matrix and r\ is a vanishingly small positive number corresponding 
to an adiabatic switching on of the field from the time t = — 00. The external field (||) induces a 
difference 5p between the neutron and proton density distributions which exhibits the same space 
time dependence 

<t\5 P \t >=ae- l(i - r e- i(uJ+ir ' )t . (2) 

The response function (or strictly speaking the retarded response function) is determined by the 
corresponding polarizability, i.e., the ratio of the density change to the field strength (which depends 
only on uj and q ): 

n(w,q) = a/e (3) 

In the case of the Skyrme effective force the retarded response function is found to be of the following 
form : 

m \ = n (u,q) 

lW,qj l-V UoM-V 1 I 1 (u,q)-V2l2^,q)' { ' 

In this formula Vq, V%, V% are related to the parameters to, t\, t2,t$, xq of the Skyrme interaction [||, [?], |8| 
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via the following formulae 



V 1 = -fa-t 1 ), (5) 
8 1 

^2 = r^(t2-5ii) , 
lb 

where po is the equilibrium density of nuclear matter and where To is the corresponding value of the 
kinetic energy density 

r(r) = £|V^(r)| 2 . (6) 

i 

In Eq. (f|) IIo is the unperturbed response function (often referred to as the Lindhard function ||): 

n (w,q) = 775— T3 / d a &^— — — (7) 

(27r) d J u + it] - e(k) + e(k + q) 

where e(k) is the single particle energy and /(k) = 1/(1 + exp/3(e(k) — /j,)) the associated occupation 
number, (5 = 1/T being the inverse temperature and \i the chemical potential. We have used the 
standard units fi = c = 1. The functions 1\ and I2 appearing in Eq. (Q) are the generalized Lindhard 
functions defined by 

T< \ 1 (aH^ /(k + q)-/(k) 

^i(^)Q) = 77: rjj / d r k.q ■— 8 

(2ttY J uj + z?7 - e(k) + e(k + q) 

w , 1 f , 3 , , 2 /(k + q)-/(k) 

I 2 (w,q) = / d J /c fc : — r-, (9) 

v M/ (2vr) 3 7 w + zry - e(k) + e(k + q) w 

Analytical expressions for the real and imaginary parts of Ho(w,q) can be found in reference §. For 

instance the imaginary part is given by 

m *2 ]_ + e /3(^-_B_) 

9m II (u,, q) = -_ log l + e ^_ E+) (10) 



where m* is the effective mass in nuclear matter calculated with the Skyrme force §, while 

m* q 2 
2m* 

The real part can be expressed as: 



E± = —2^±7^Y- (11) 



„ „ /" m*/c f fc / , ,m*u; o . , ,m*uj n A] , m . . . 

,ien »=-y^{- i+ ^(* ( i5- + i ) -^-iv} mkx> <i2) 

with: 

0(x) = (l-x 2 )log . (13) 
3; + 1 
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From the previous equations the value of I\ can subsequently be calculated from the identity 



m*(w + e(q))n + 2Ii(w,q) =0 . (14) 

In the case of zero temperature, expressions for I 2 can be found in Ref. ||. For non zero temperature 
the calculation of the integral (^) is similar to that of Ho. The result is that the real part of I 2 is an 
average of the zero temperature values calculated for the same energy uj and momentum q, but with 
various Fermi momenta k F distributed with a weight factor which is just the derivative of the Fermi 
occupation number. Explicitely on has the following formula: 

I 2 {uj,q,T) = -J 3fie I 2 {uj,q,T = 0,k F = k) df(k,T). (15) 

The imaginary part of I 2 is given by 

Qml 2 (uj,q) = - 27rq p2 X > ( 16 ) 

with 

1 4. p P(ji-E-) 

X = (3E + log l + e ^_ E+) + Li 2 (l + e^' E ^) - Li 2 (l + ^~ E -)), (17) 



where Li 2 the Euler dilogarithmic function |ll| and E± has been defined in (11). 

The imaginary part of II (c.f., Eq. (|j)) is related to the value of the strength per unit volume S(lu) 
for the operator T3 exp(iq.r) by 



S(u) = 3mll(uj,q). (18) 

7T 



It satisfies the following sum rule 



00 ^,2 







S(u)ucL; = 7 f 1 ;po(l + K) (19) 
2m* 



where k is the enhancement factor due to the momentum dependent terms of the Skyrme interaction: 

2mkn 



ih + h) (20) 
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The strength distribution as a function of the frequency uj is plotted for the three Skyrme forces SV, 
SI and SII in figures 1, 2 and 3 respectively, for various values of the temperature. In these figures 
the value of the momentum transfer was chosen to correspond to the dipole mode in lead-208 in the 
Steinwedel- Jenssen model. In this model neutrons oscillate against protons inside a sphere of radius 
R (the nuclear radius) according to the formula ||] 

< r|<^o[r >= esin(q.r) sin(u;£), (21) 
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where 

" = m (22) 

The radius of lead-208 is i?=6.7 fm. This value leads to g=0.23 fm . 

It can be noted in Fig.l that in the case of the Skyrme force SV, the dipole strength is nearly 
independent of temperature. We found this to be the case for all Skyrme forces having small value 
of the effective mass, such SV (m*=0.36, in units of the bare nucleon mass whenever a numerical 
value is given) and SIV (m*=0.47) @. These forces also give an energy of the dipole resonance in 
lead-208 much higher (about 25- 30 MeV) than the observed value (15 MeV). In contrast effective 
forces with an effective mass of the order of unity such as SI (m*=0.91) and SVI (m*=0.95) ||, give 
a sharper resonance at zero temperature, a good agreement with the observed position and also a 
stronger temperature dependence. This result may be related to the saturation of the multiplicity 
of photons emitted by the giant dipole resonance in hot nuclei when their temperature exceeds T=3 
MeV observed in some experiments | jPIJ| . 

Concerning the collectivity of the resonance one can see from Fig. 2 that the SI Skyrme force 
gives a more concentrated coherent behavior. Taking forces with smaller value of effective mass, for 
example SII (m* = 0.58) in Fig. 3, the interaction ceases to be important and the response function 
has pratically the form of the bare response, i.e., the imaginary part of the unperturbed Lindhard 
function. With a smaller effective mass ( as in SIV and SV), the collectivity is distributed over an 
even larger range of frequencies. 

It is worthwhile emphasizing that all Skyrme interactions considered here provide a good descrip- 
tion of nuclear ground state properties such as root mean square radii and binding energies. Obviously 
such a good description is not sufficient to garantee at the same time an accurate description of the 
isovector dipole response function: even with a good fit to the symmetry energy coefficient there is a 
large dispersion in the results for the nuclear matter isovector response function. 

Our main conclusion is thus that, even at the level of the random phase approximation, uncertain- 
ties in the predictions of microscopic calculations do remain for the dipole strength and its temperature 
dependence. We find that the effective mass mainly determines the energy of the resonance, its col- 
lectivity and its evolution with temperature. However other effects are known to be important in 
the description of the observed width such as the coupling to two particle- two hole configurations 



\12\. A larger number of experimental data on the temperature dependence of the strength would be 
highly desirable to resolve the remaining uncertainties. Simultaneously further microscopic theoretical 
studies of the effects not included in RPA would be necessary. 
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